Sagnac interference in Carbon nanotubes 
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The Sagnac interference mode arises when two interfering counterpropogating beams traverse a 
loop, but with their velocities detuned by a small amount 2u, with Vh/l = vf iu. In this paper we 
perform a perturbative non-equilibrium calculation of Sagnac interference in single channel wires as 
well as armchair nanotube loops. We study the dependence of the Sagnac conductance oscillations 
on temperature and interactions. We find that the Sagnac interference is not destroyed by strong 
interactions, but becomes weakly dependent on the velocity detuning u. In armchairs nanotubes with 
typical interaction strength, 0.25 < g < 0.5, we find that the necessary temperature for observing 
the interference effect, Tsag is also only weakly dependent on the interaction, and is enhanced by a 
factor of 8 relative to the temperature necessary for observing Fabry-Perot interference in the same 
system, Tfp- 
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I. INTRODUCTION 



One of the most tantalizing effects predicted by quan- 
tum mechanics is the appearance of interference fringes 
when two matter beams come together. These fringes 
provided the ultimate testimony to the pertinence of 
quantum mechanics and the Schrodinger equation. Inter- 
feromctry of light is employed in many precision measure- 
ment devices. The Mach-Zchnder interferometer pro- 
duces interference between two beams traversing two dis- 
tinct paths, one of which passes through a test cham- 
ber containing, for instance, a dilute gas (see Figfl]); this 
setup was originally used to measure the refraction index 
of the gas in the chamber. Fabry-Perot interferometer 
recombines a scries of beams, where the n'th beam tra- 
verses the optical path between two mirrors or through 
a loop n times. The narrowness of the resulting inter- 
ference peaks allows a precise measurement of a light 
beam's wave length, and is commonly used to measure 
the Zeeman splitting of an atom in a magnetic field. The 
most sensitive of all interference constructs, however, is 
the Sagnac interference^. In this setup, a light beam is 
split into two beams, which traverse the interferometer's 
loop both clock wise and counter clock wise, before be- 
ing recombined. In this case, the interference fringes arise 
due to an absolute rotation, and provide the most accu- 
rate measure of the angular velocity of the device. This 
was used by Michclson to measure the absolute rotation 
of the Earth. More recently, the Sagnac interference ef- 
fect was cleverly used to measure time reversal symmetry 
breaking in superconductors 2 -. 

Quantum mechanics opened the way for matter-wave 
interferometry. Electron interferometry is a powerful 
probe of interaction effects on low-energy phases of quan- 
tum matter, as demonstrated by numerous examples. 
Mach-Zchnder interferometers reveal Aharonov-Bohm 
oscillations and quantum hall effect edge channels^£&£, 
and can probe exotic fractional quantum Hall state a 8 ' 9 ' 10 . 
Similarly, two-path Mach-Zchnder interferometers can 
probe correlated states of quantum dot o 11 ' 12 . Of par- 
ticular interest to us are metallic carbon nanotubes. The 
Luttinger liquid behavior in these system a 13 ' 14 i 15 was 



partially verified through the observation of Fabry-Perot 
interference in finite sections of the nano-tubei^. The 
Fabry-Perot interference should, in principle, allow the 
observation of spin-charge separation and determination 
of the interaction parameters of the Luttinger liquid^. 
But the similar energy scales of the spin and charge 
modes' interference patterns has made such experimental 
observation challenging. 

The most sensitive interferometer of all, however, the 
Sagnac interferometer, has not been seriously explored 
yet in the context of interacting electronic systems. In 
Ref. [H we proposed that this interference naturally 
occurs in metallic armchair nanotube loops (Fig. [2]). 
Instead of rotation, the Sagnac interference arises due 
to the band velocity difference between right- and left- 
moving electrons about each Dirac node. This velocity 
difference is present whenever the electronic Fermi sur- 
face is tuned away from the Dirac points at half-filling, as 
shown in Fig. [3ja). The operating principle of the elec- 
tronic Sagnac effect has the same origin as the universal 
conductance fluctuations, and weak- localization effects 
in disordered two-dimensional electron gase a 19 i 20 ' 21 . In 
nanotubes, it can also appear due to band-scattering in 
a pair of impurities^ 2 -. 

Because the Sagnac effect involves electrons traversing 
the same path in two different directions, rather than re- 
peating the same path as in Fabry-Perot interference, the 
phase accumulation is extremely small. Therefore Sagnac 
interference exhibits large-period conductance fluctua- 
tions as a function of gate- and source-drain voltages, and 
is expected to persist to high temperatures in comparison 
to Fabry-Perot interference, which is more sensitive to 
thermal dephasing. This interference mode should thus 
be able to reveal much more precise information about 
the unique state of interacting electrons in thin quantum 
wires. 

Our goal in this manuscript is to thoroughly explore 
the range and robustness of the Sagnac interference 
mode, concentrating on armchair Carbon nanotubes. 
The questions we will ask concern the amplitude of this 
interference mode as a function of the temperature, gate 
and source-drain voltage, and Luttinger parameter of the 
nanotube. 
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a. Mach-Zehnder 




FIG. 1: (a) In a Mach-Zehnder interferometer the input beam 
is split into two beams which traverse independent paths be- 
fore being recombined. (b) In a Fabry-Perot interferometer a 
beam is split into a deflected ray, which is recombined at the 
output with a ray that traverses a loop, (c) The Sagnac in- 
terferometer splits the beam into a two beams which traverse 
the loop in two opposite orientations, and get recombined at 
the output. This allows a very sensitive measurement of the 
angular velocity of the interferometer, as it results in a differ- 
ent relative speed in the clockwise and couterclockwise rays. 
Clear rectangles represent beam-splitters, and patterned rect- 
angles represent mirrors. 




FIG. 2: Schematic of a nanotube loop reported in Ref. Il8l . 
This geometry allows electrons to tunnel from the point X 
on the loop to a distant point X' on the other end of the 
loop, and vice versa. We refer to this process as cross-loop 
tunneling. An electron entering from the left can traverse 
the loop moving right with velocity vr, without scattering, or 
tunnel from X to X' and traverse the loop moving left with 
velocity vl- 
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FIG. 3: (a) The energy spectrum of an armchair nanotube. 
When the chemical potential is tuned away from the degener- 
acy points by a gate voltage, the left and right movers in each 
node will have different velocities, which leads to the Sagnac 
interference in the loop geometry, (b) The scatterings we con- 
sider in this paper tunnel, say, a right moving electron from a 
given node at point X to a left moving electron, of the same 
node, at point X' , and vice versa, (c) Sagnac interference can 
also arise without the loop geometry through inter-node tun- 
neling, since right movers at node 2 have the same velocity as 
left movers at node 1. 



II. SAGNAC INTERFERENCE IN A SINGLE 
CHANNEL 



The paper is organized as follows. In Section [TH as a 
warm-up, we analyze the simpler case of Sagnac inter- 
ference in a single channel of right- and left-moving elec- 
trons. In III Al we introduce the model of a single channel 
with a linearized spectrum, and the cross-loop tunneling 
which will give rise to the Sagnac interference. In III Bl 
and III Cl we set up the non-equilibrium perturbative cal- 
culation of the conductance in the presence of cross-loop 
tunneling, and in III Dl and III El we analyze the behavior 
of the oscillating conductance as a function of gate and 
bias voltages and temperature. In Section Mil we repeat 
the above steps for the physically relevant case of Car- 
bon nanotubes, including spin and node degeneracies in 
the calculation, and remark on the similarities and differ- 
ences from the single channel case. Finally we conclude 
with a discussion of the experimental implications of our 
calculations. 



As discussed in the introduction, the Sagnac interfer- 
ence in the loop geometry is due the the asymmetry be- 
tween the velocities of the left and right moving electrons. 
To demonstrate this in the simplest form, we first study 
in this section a single channel with a single type of left 
and right movers. In a carbon nanotube, there will be 
four such channels due to spin and node degeneracies. 



A. The Model 

We start with a single one dimensional channel of elec- 
trons and a linearized spectrum, with different left and 
right mover velocities, and a density-density interaction. 
The Hamiltonian density for this system is: 
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+ a (ii>y> R + i>y L 



(i) 



where the operator ^ R i L creates a right /left moving elec- 
tron, with the velocity: 



vr/l = v f ±u. 



(2) 



The scattering we are interested in is the one which takes 
a right moving electron at one side of the loop, point 
X is Figure [H and scatters it to a left moving electron 
at the other side of the loop, point X 1 , and vice versa. 
This process has been dubbed Cross-Loop scattering in 
Ref. [lg]. The same effect can also be obtained without 
the loop geometry by inter-node tunneling^, since right 
movers at node 1 move with the same velocity as left 
movers at node 2. This inter- node tunneling is shown in 
Figure [3|c). If we choose our coordinate along the loop 
such that the point X corresponds to x = and the point 
X' corresponds to x = L, then this scattering process is 
described by the Hamiltonian: 



H ba =T^ R {Q)ML) + h.c. 
+T 2 ^ L (0)^ R (L) + h.c. 



(3) 



In the presence of the quartic density-density interac- 
tions in the Hamiltonian, Eq. (JT|) , it is useful to use the 
standard bosonization procedure , since the Hamiltonian 
is quadratic in terms of the bosonic fields. The electron 
fields are bosonized as follows: 



il>R/L 



J(<t>±6) 



(4) 



where 9 and <f> are bosonic fields that satisfy the com- 
mutation relations [0(x), 4>(x')] = i(jr /2)sgn(x~ x'); also, 
the total density and the current density are given by 
iV6> = pr + p L , and ±V</> = p R - p L , respectively. The 
Hamiltonian in terms of the bosonic fields becomes: 



H\ c h = — / dx 
2tt 



-(we) 2 

9 



■g(V<j>y + 2-V(9V0 
v 



(5) 
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where q - 
parameter and v = Vp/g. 



-1/2 



is the Luttinger interaction 
This is the familiar Hamilto- 



nian of a ID interacting electron system, with the ad- 
dition of the u term which gives left and right moving 
particles different velocities. Indeed, this Hamiltonian 
can be easily diagonalized and the left and right veloci- 
ties turn out to be for a general value of the interaction 
parameter g: 



Vp 

Vr/l = V ±u = — ± u. 

9 



(6) 



Our goal is to calculate the effects of the Sagnac in- 
terference as seen in the conductance as a function of 



the applied bias and gate voltages, and as a function of 
temperature. Due to the applied voltages the system 
is not in equilibrium, and we must turn to the Keldysh 
non-equilibrium formalis m 23 i 24 . Below we carry out this 
analysis first for the simplified electron gas with the scat- 
tering Hamiltonian Hb s , Eq. ©, as a perturbation. 



B. Non-Equilibrium correlation functions and 
conductance 

The response of the loop to a bias source-drain volt- 
age can be analyzed using the non-equilibrium Keldysh 
formalism. Following RefUTI. we assume that in the dis- 
tant past, before turning on the backscattering, the left 
and right moving electrons separately had well defined 
thermal distributions set by separate chemical potentials. 
The density matrix corresponding to this initial distribu- 
tion at temperature T = 1//3 is: 



Pv 



1 



(7) 



with Zy = Tr[e 13 Hv ] and the Hamiltonian which takes 
into account the applied voltages is: 



Hi 



■h 



H 



e— (Nr 

V sd 

left - e-— 



N L )-aeV g (N R + N L ) 



dx 



V0 



aeV n I dx 



V0 



(8) 



The gate voltage, V g , simply couples to the total charge 
density, with a being a geometrical factor of the sys- 
tem, while the source-drain voltage, V s d, induces the im- 
balance in the chemical potentials of the left and right 
movers. 

As explained in Ref. [I?], both V s d and V g can be elim- 
inated from the unperturbed action by an appropriate 
unitary transformation, which is equivalent to shifting 
the bosonic fields by a function of space and time; this is 
easy to see if one writes down the Lagrangian including 
the voltages^. The equivalent shifts for the case at hand 



ag 2 eV g 



1 



hvp 1 — g 2 u 2 /v 2 
ag 2 eV g u/vp 
hvp 1 — g 2 u 2 /iij 



eVsd 

2h 1 



(9) 



These shifts remove the voltages from the Hamiltonian 
Hy and therefore all the correlations to appear in the 
calculation will be equilibrium correlation functions with 
respect to H\ c h ■ The dependance on the applied voltages 
now appears in the scattering Hamiltonian, Hb s , due to 
the shifted bosonic fields. 

Let us now focus our attention at the charge current, 
which in the bosonic language is I = (e/n)dtO. After 
performing the unitary transformation described above 
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FIG. 4: The Keldysh contour used in the non-equilibrium 
calculation. The Keldysh time ordering operator Tk orders 
operators along the contour, so fields on the (+) branch is 
always at an earlier time than fields on the (-) branch. 



space, these correlations arc functions of x — x' and t — t', 
for example: 



C?+(x,t) = C°*(x,t',0,0) 

log w fl smh 

V V P V R 



log ( vl sinh 



(vLt + x)n 
(3vl 



(13) 



we can write the formal expression for the expectation 
value of the current in the usual interaction picture^: 

(/) = to+^—Tr (er^f K {l K (x,t)e-*fc *^.<f)}) 

(10) 

Tk is the time ordering operator along the Keldysh 
contour shown in Fig|4l and lK{x,t) is the symmetrized 
current operator with respect to the two branches of the 
contour. The current Iq = e 2 V s d/h is the ideal current 
that would flow in the absence of backscattering in a com- 
pletely transmitting channel, and it explicitly appears 
due to the shift of the 9 field. The Hamiltonian H' bs de- 
notes the scattering Hamiltonian Hb s with the properly 
shifted bosonic fields. The expression for the current can 
be expanded in powers of H' bs , and all the correlation 
functions to appear in this expansion are equilibrium cor- 
relation functions at temperature 1/(3. If we denote by 
9 + and 6~ the fields on the forward branch and back- 
ward branch of the Keldysh contour respectively, then 
time ordering along the contour means that 9 + 9 + corre- 
lations have the usual time ordering, 9~9~ are anti-time 
ordered, and 9 + (t) is always earlier in time that 6~{t'). 
The same applies for all the fields. 

It is useful to apply a Keldysh rotation to the fields, 
9 ± = 0±|#, and similarly for </>. The correlation function 

(TK@(t)6(t')} vanishes by construction, and we define : 
C 9 (x, t; x', t') =(T K 6(x, t)6{x',t')) 
=±{{9(x,t)J(x',t')}) 

R e {x, t; x',t') ={T K 6(x, t)9(x', t')) 

= -ie(t-t')([§(x,t)J(x',t')}) (11) 

and similarly for the 4>j fields, and for the mixed corre- 
lations: 

C e<t> {x, t; x', t') =(T K 6(x, t)4>{x' ,t')) 

=±({9(x,t),4>(x',t')}) 

R^ix, t; x', t') =(T K 6(x, t)4>(x', t')) 

= -ie(t-t')([6(x,t),4>(x',l/)}) (12) 

where operators with a hat are simply the time dependent 
operators with no time ordering. As explained above, 
these correlation function are to be evaluated in equilib- 
rium, and therefore are easily explicitly calculated (Ap- 
pendix [B]) . Due to translational invariance in time and 



C. Perturbation Theory 

The Sagnac interference fringes occur already with 
weak bakescattering at the base of the loop, and can 
be deduced from a perturbation analysis of the tunnel- 
ing Hamiltonian, Eq. ([3|). As outlined above, to calcu- 
late the current, I\ c h = (^dt9), we absorb the gate and 
bias voltages, V g and V s d respectively, in the shifts in 
EqlQl which allow us to move the voltages from the un- 
perturbed Hamiltonian H\ c h to the backscattering per- 
turbation, Hi, s . Then, we expand the formal expression 
we found for the current using the Keldysh technique, 
Eq. [lOl in powers of Hb Sl and use Wick's theorem to 
evaluate the resulting contributions. 

To lowest nontrivial order, which is second order in 
Hf, s , we obtain after a lengthy calculation: 

e 2 V sd 

Ilch = ; h I co + hnco (14) 

h 

The first term is simply the current that would flow 
through the system in the absence of backscattering. The 
coherent current, I col oscillates with the gate voltages V g , 
and is given by: 

/ 2ug 2 La \ 
I co = cri r 2 cos \-^-^—^ Vg j x 

J dt sin(^t) e - c "(i,*) si n (i? co (L ; t )) (15) 

where c is a constant of order unity, and we assume that 
Tj are real for simplicity. The incoherent current, Ii nco , 
is independent of the gate voltage, and is given by: 

c f r ?E / dt^{^t)er c ^ L ^^{R mc AL^ 

+ c (ri ^ T 2 ,L — > —L) (16) 

The functions C co , Cf nco , R co and Rinco are compli- 
cated combinations of the correlation functions defined 
in section III Bl and are given explicitly in Appendix [B] 
These functions do not simplify, partly due to the fact 
that the correlation functions in this problem are not 
symmetric under x — > — x since left and right movers 
have different velocities. 
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FIG. 5: Differential conductance oscillations of a single chan- 
nel of fermions, as a function of bias voltage V s d> f° r velocity 
detuning u/vf = 0.1 and interaction strength g = 0.5. The 
beating is due to the only two voltage oscillation frequencies in 
the problem, f2i = , , eL , and 0,2 = * r eL , , where v = vf/q- 

1 1 h{v+u) n(v — u) ' J / 3 

The voltage is in units of HvF/eL. The shorter voltage os- 
cillation periods is AV sd = 2tt(^^±^.)~ 1 « 12.5hv F /eL, 

and the large oscillation period is AV a d = 27r ( ni ~° 2 ) 1 
250 /eL. 

D. Voltage dependence of the single-mode Sagnac 
interference 

The voltage current characteristics given in Eqs. ([HI 
- [16]) can be evaluated numerically to obtain the voltage 
and temperature dependence of the single-mode Sagnac 
interference. The period of the interference as a function 
of the gate voltage (I co ) are easily observed to be (for 
small u/vf)'- 

Av s agnac ^vrnVvr^ 

u ag z L u y 

where ^V pp is the period in gate voltage for Fabry-Perot 
interference. Fabry-Perot interference occurs whenever 
part of the wave's trajectory can be repeated. Since the 
Sagnac interference involves traversing the same path in 
two different directions, the phase difference accumulated 
in the process is much smaller than the difference in- 
curred by repeating part of the path, and therefore the 
period of the Sagnac interference is much larger than the 
period of the Fabry-Perot interference. Such large period 
oscillations have been experimentally observed in Carbon 
nanotubes, in the loop geometry, as reported in Reffl8l. 
in addition to the shorter period Fabry-Perot oscillations. 

For a given gate voltage, both the coherent and inco- 
herent parts of the current oscillate with the bias voltage 
V s d- This oscillation is due to the fact that in the pres- 
ence of bias voltage, the Fermi energy of the left- and 
right- moving electrons are different by V sc i, and hence 
their Fermi wavevectors are different also and they would 
acquire different phases traversing the loop. This os- 
cillation will be present even for no velocity detuning, 



u = 0. When the detuning is finite, u ^ 0, the differ- 
ential conductance G\ c h = d!ich/dV s d will show a beat- 
ing pattern due to the two different left and right mov- 
ing excitation velocities. Here we are only considering 
the Sagnac oscillations arising from the cross-loop tun- 
neling, Eq. © . Figure [5] shows the oscillations of the 
differential conductance at a fixed gate voltage. For non- 
interacting electrons, the beating pattern corresponds to 
the addition of two harmonics with two different fre- 
quencies in voltage, sin^^V^) and sm^^V^), with 
Qr,/l = hvR /L an< ^ v R/l = v f i u. The beating pat- 
tern will then display fast oscillations with voltage period 
AV^ ast = 2tt ^ [ga+fj^l j ^ anc j s i ow voltage oscillations 

with period AV s s j ow = 2ir For interacting 

fermions, g 1, the oscillations will not be simple har- 
monic oscillations. Still, the periods will be evident and 
will have the same functional form, in terms oivn/L, as 
the frequencies in the non-interacting case. The periods 
AV/ d ast and AV^ ow do depend on g through velocities 
v L and vr, v R/L = ± u (Eq. [6]). 

These oscillation, generally, lie atop a powcrlaw be- 
havior of the differential conductance as a function of 
V s d, as expected from the known behavior of the conduc- 
tance in the presence of impurity backscattering2£. For 
backscattering from an impurity in a Luttingcr liquid, 
the backscattered current, for low temperature, behaves 
as / cx V^ 2 d 9_1 . We chose to plot the Sagnac oscillations 
as a function of V s d (Fig. [5]) for the interaction parame- 
ter g = 0.5 since for that value the corresponding power 
law would be / oc V^, and the contribution of such a 
powerlaw to the differential conductance would vanish, 
making the oscillation atop this powerlaw more visible. 

E. Temperature dependence of the single-mode 
Sagnac interference 

Next we consider the temperature dependence of the 
gate- voltage driven oscillations in the coherent part of the 
current. As argued in Rcf. [3, the large period Sagnac 
oscillations are expected to be observed at much higher 
temperature than the shorter period Fabry-Perot oscil- 
lations. This difference in temperature behavior can be 
easily understood by examining the phase giving rise to 
the interference in both cases. In the Fabry-Perot case for 
a loop, the lowest order interference is between a beam 
of electrons which is not scattered, and a beam of elec- 
trons which, due to scattering at the base of the loop, 
does a roundtrip between the the two scattering points. 
The phase difference between these two beams at energy 
E is Aippp = kpjL = -J-L-^. Finite temperature effec- 
tively causes uncertainty of order T in the energy E, and 
the interference pattern will be washed out when the un- 
certainty of the the phase A(/)fp is of order 2ir, which 
happens at a temperature Tpp = ^j^vr. 

In the Sagnac case, the interference is between a beam 
that traverses the loop moving left and one which tra- 
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verses the loop moving right. The phase difference 
between these two beams at energy E is A<Psag = 

kj L — kr>L = ( — -]L|, and this interference will 

-l 



be washed out at temperature Tsag = — 

For non-interacting electrons the right and left moving 
velocities are = Vp ± u. Thus to lowest order in 

u/vp, the highest temperatures for observing interference 
according to the argument above are: 



Ti 



FP 



TsAG 



ttHvf Vp 



FP 



VF 

u 



(18) 



For non-interacting electrons, we expect the Sagnac inter- 
ference to survive to a temperature higher by a factor of 
vpju than the corresponding Fabry-Perot temperature. 
We will show through explicit calculation that this is in- 
deed true for the non-interacting case. For interacting 
electrons, we will see that Tsag will still be considerably 
larger than Tpp, but their ratio is less than the dramatic 
vpju ratio. 

To explore the Sagnac temperature range, we evaluate 
the amplitude of the coherent oscillations (the oscilla- 
tions in V^) as a function of temperature, for different 
interaction parameters g and different ratios of u/vp. 
For non-interacting electrons, g = 1, we find that the 
Sagnac oscillations indeed survive up to a high temper- 
ature, which is a factor of vp/u higher than the cor- 
responding Fabry-Perot oscillations. Figure [5] plots the 
oscillation amplitude as a function of temperature, nor- 
malized by its zero-temperature value, and for different 
values of u/vp. The functional dependence on tempera- 
ture is given approximately by: 



G co {T) 
G co (T = 0) 



2nk B LT 

hvp 



Vp 



1 



sinh^Tr/cBT-p^-^) 



(19) 

This result is similar to the exact form of the temperature 
dependence of the Fabry-Perot interference amplitude^, 
with the only difference being the factor oiu/vp. There- 
fore, the Sagnac oscillations of non-interacting electrons 
indeed survive up to a temperatures which are a factor 
of vf ju larger than the Fabry-Perot oscillations. 

For interacting electrons, j ^ 1. the Sagnac interfer- 
ence still survives up to temperatures significantly higher 
than the corresponding Fabry-Perot temperature scales, 
but the enhancement is suppressed compared to that of 
non-interacting electrons. Figure [7] shows the Sagnac 
temperature scale T* vs. u/vp for three different values 
of the interaction parameter g, where we define T* to be 
the temperature at which the amplitude of the oscilla- 
tions reaches e _1 of its amplitude at zero temperature. 
For non-interacting electrons T* is strongly dependent 
on the ratio ujvp as discussed above. For the interaction 
parameter values g = 0.5 and g = 0.25 (Dashed lines), 
the temperature T* is only weakly dependent on the ra- 
tio u/vp. As an example for the resulting enhancement 
of the Sagnac compared to the Fabry Perot interference, 
consider g = 0.25, where the T* temperature scale for the 
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FIG. 6: Coherent Sagnac oscillation amplitude vs. tempera- 
ture for different values of (vf/u), for noninteracting elec- 
trons (g — 1). The slowest decaying plot corresponds to 
vf/u — 100, and the fastest decaying plot corresponds to 
vf/u=10. Temperature is given in units of Flvf /ksL. 
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FIG. 7: T* vs. u/vf, where T* is the temperature at which 
the coherent differential conductance (the part of the con- 
ductance which oscillates with gate voltage) reaches e _1 of 
its zero temperature value. For a non-interacting system, 
g = 1, the single channel case gives the same temperature 
dependence as the case with spin and node degeneracies, 
T* oc vf/u. The single channel temperature dependence is 
given for g — 0.5 (squares, dashed), and g = 0.25 (diamonds, 
dashed). The Carbon nanotube temperature dependence is 
given for g = 0.5 (triangles), and g = 0.25 (inverted triangles). 
Temperature is given in units of ft«i?/fcsL. For reference, the 
T* corresponding to the g — 1 Fabry-Perot oscillations is also 
plotted. 



Sagnac oscillations is roughly 1.6hvp /kp,L, a factor of 4 
enhancement over T* of the non-interacting Fabry-Perot 
oscillations which is 0A2hvF/kE,L, despite the suppres- 
sion of the Sagnac T* due to interactions. As can be 
seen in the figure, for g = 0.5 the enhancement is about 
7. While it is difficult to extract the analytic dependence 
of the temperature on the interaction parameter, one can 
repeat our calculation for any value of g. 
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III. INTERFERENCE IN NANOTUBES 

Equipped with our understanding of the single-channel 
Sagnac interference, we can now consider the likely phys- 
ical system where it may be observed: a metallic Carbon 
nanotubc with four different Dirac nodes. We now add 
the spin and node degeneracies of a Carbon nanotubc, 
and examine their effect on the Sagnac interference pat- 
tern voltage and temperature dependence. 



Next we bosonizc the electron field operators in the 
nanotube. The slowly oscillating parts can be written 
as: 



1pR/La 



(23) 



where 9 aa and (f> aa are bosonic fields that sat- 
isfy the commutation relations [9 aa {x), 4> a 'a' = 
i(n/2)5 a ,a'8a.a'sgn(x — x'). The Hamiltonian in terms 
of the bosonic fields is 18 : 



A. The Model 

The energy spectrum of a Carbon nanotubc is shown 
in figure [3] a. This spectrum is usually linearized around 
the Fermi surface, which yields four chiral modes, two left 
moving and two right moving (not including spin), with 
linear dispersion. These modes can be bosonized and 
treated within the Luttinger Liquid theory framework, as 
we have done in the single channel case in the previous 
sections. All these modes are usually assumed to have the 
same velocity, the Fermi velocity vf- For the purposes 
of this paper, it is important to notice that when the 
Fermi surface is away from the degeneracy points where 
the upper and lower bands meet, linearizing the spectrum 
actually gives two different velocities which we shall note 
v± = vp±u. The linearized Hamiltonian density is, then: 

2 

Ti-ich = «E E { V Ra^ Raa d x ^ R aa ~ V La l\)\ aa d x l\) Rac 
a=l cr=T4 

2 



A 



E E (^L>flaa + ^La^Lac 
a=l<r=T,J. 



(20) 



where ipR/Laa stands for a right /left moving electron at 
node a with spin cr, and we added a total charge den- 
sity interaction term. The velocities that appear in the 
Hamiltonian arc: 



VR/Lltr = VF ± U = V± 
VR/L2a = VF T U = V T . 



(21) 



Thus for u > 0, at node 1 right movers are faster than left 
movers, while at node 2 the opposite is true. Now, the 
nonlincarity of the clctronic spectrum in a Carbon nan- 
otube needs to be taken into account when considering 
the velocity difference, u\ it depends on the detuning of 
the chemical potential away from the degeneracy points. 

The scattering process we are interested in is very sim- 
ilar to the one we had in the single channel case. We need 
to consider a term that scatters a right mover at one end 
of the loop to a left mover at the other end of the loop, 
conserving spin and node quantum numbers, 



E [ r iV>ka(0)Vw(£) 



h. 



a,a=l,2 



+r 2 V>L<r(o)VW£) + h. 



(22) 



H 



Ach 



2-K 



,a=1.2 L 



+ (-l) o+1 2£V<fc I(T V0 a(7 + / dx\ 



E £V<£ 

er,a=l. 2 



(24) 

If the velocities of all branches of the spectrum were 
equal, i.e. u = 0, then the Hamiltonian Hb would be 
diagonalized by the spin and node symmetric and an- 
tisymmetric combinations of the 8's and <fi's~. By di- 
agonalizing we mean a linear mapping of the <p and 9 
fields such that the Hamiltonian takes the form of four 
independent channels, each resembling of H lc h, Eq. |5]). 
When u ^ 0, there still exists a local transformation 

&aa = Ej=i..4 { AJ aa^j + B la^j) tnat diagonalizes the 
Hamiltonian, but it is a more complicated combination of 
the fields that depends on u and A, and mixes the theta 
and <p fields, which makes the conductance calculation 
quite cumbersome. While the details of this transforma- 
tion are given in appendix El the diagonal Hamiltonian 
is: 



H ^ h = E § I dx 



j=3,4 



E ^ / 



J=l,2 



2tt 



-(v^) 2 + . 9j (v^) 2 

9j 



+ (-l) :i+1 2—V(j) j V9 J 

Dp 



(25) 



The fields 9i/ 2 and 4>i/2 are t ne spin antisymmetric com- 



binations of ( 



1/2, 



and 



, l/2< 



respectively. Since the inter- 



action term in Eq. (|20[) involves only the spin symmetric 
combinations, the spin antisymmetric combinations are 
untouched and still have the left and right moving ve- 
locities as in Eq. ([2T|) . On the other hand, the fields 6*3/4 
and </>3/4 are not simply the remaining symmetric com- 
bination and mix the remaining 6>'s and </>'s. These fields 
have the same left and right moving velocity, which is: 



V3/4 



v F 

vn 1 




where g = ( 1 + J is the Luttinger parameter^ 

Fortunately, for the region of parameters which is 
of interest, namely strong interactions, g < 0.5 and 



(26) 



ujvp < 0.1, the exact change of basis required to di- 
agonalize the spin symmetric part of the Hamiltonian 
is very close to the usual node symmetric/antisymetric 
basis. This can be explicitly seen, for example, from the 
velocities of these modes. For this entire range of parame- 
ters, the velocities of the diagonal fields, given by Eq. (|2l))) . 
are at most 1% different from the values wc expect for the 
left-right symmetric system, which arc vp/g and vf- Due 
to the strong interactions in this spin symmetric sector, 
the velocity asymmetry is unimportant, and it is for this 
reason that we choose to still use the node symmetric- 
antisymmetric basis and treat these fields as the diag- 
onal ones. In Appendix [A] we elaborate on and justify 
this approximation. Note that the velocity asymmetry is 
still apparent in the non-interacting spin antisymmetric 
modes labeled by j = 1 and j = 2 in Eq. (|2"5]l. 



B. Perturbation Theory 

Using the diagonal form of the Hamiltonian with the 
above approximation, wc proceed to calculate the cur- 
rent, I = (e/n)(Ylao- dt^aa), as in Section |TTJ The ap- 
plied voltages now couple to the total density and total 
number of left movers and right movers: 



Gtube {arb) 



(a) 



H v = 



v sd 



(Nr -N l )- aeV g (N R + N L ) 



(27) 



The external voltages can be removed from the Hamilto- 
nian by the appropriate shift of the bosonic fields: 



o.oi 




50 




0.025 
0.020 
0.015 
0.010 
0.005 



Fourier[G tU 6e] 



(b) 




1 ftGfigr) 

FIG. 8: (a) Differential conductance oscillations for a nan- 
otube, i.e. including both spins and both nodes in the spec- 
trum, for velocity detuning u/vf = 0.1, and interaction 
strength g — 0.5. In the nanotube case, the beating is due to 
the four voltage frequencies in the problem, fli = S^, where 
v i/2 = vf i u, vs ~ vf and 114 m vfIq- (b) The voltage 
Fourier transform of the oscillations in (a) clearly displays 
the four dominant frequencies , fii_4, corresponding to the 
four velocities in the problem, and encode the nanotube pa- 
rameters vf, g and u/vf- 



ag 2 eV g 



1 



hvf 



1 



g 2 u 2 /v 2 



«w + (-i) 



a ag 2 eV g 



F 

u/vf 



eV s d 

2h 



Hvf 1 — g 2 u 2 /v 



.,2 • 
r 



(28) 



Again we use the Kcldysh contour to write the formal 
expression for the current, as in Eq. (|10p . and expand 
it to lowest order in the appropriate H' bs which contains 
the voltage dependence due to the shifts of the fields. 
The approximation we made above, namely that it is the 
node symmetric/antisymmetric combination which diag- 
onalize the Hamiltonian, allows us to write the current 
in a very similar form to the single channel case: 



hch = 4- 



> sd 



(29) 



which oscillates with the gate voltage: 

2ug 2 La 



i rn =crir 2 cos 



h 2 v 2 F (l-g 2 u 2 /v 2 F ) 9 



V a X 



dt Bin(^t) e- d °°( L >V sin(i? co (L, t)) (30) 



and the third term is the incoherent current, which is 
independent of the gate voltage: 



c (ri^r 2 ,L^-L) 



eV s d^ --cv neo (L,t) 



sva{Ri nco {L,t))\ ^ 
(31) 



The first term on the right hand side of Eq. (|2"9l is the 
current that would flow in the nanotube in the absence of 
backscattering. The second term is the coherent current 



The function C co , C^ co , R co and Ri nco are related to 
the single channel correlation functions as explained in 
Appendix [5] 
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C. Temperature and Voltage Dependence in 
Carbon Nanotubes 

As in the single channel case, we find there is a coher- 
ent part of the interference current which oscillates as a 
function of the gate voltage with a large period, much 
larger than the Fabry-Perot oscillation period, as seen 
explicitly from the voltage dependence of I co . 

The differential conductance dl4 C h/dV s d, on the other 
hand, displays a beating pattern, but a more complicated 
one than in the single-channel case, since there are four 
different velocities in the problem now: v F ± m, 1)3 ~ vf 
and t>4 ~ vf /g. Figure [8] shows the differential con- 
ductance of the nanotube, dIi c h/dV s d, and its Fourier 
transform. From the Fourier analysis we see that clearly 
there are four dominant frequencies, which correspond 
to the four different velocities of the collective modes in 
the nanotube. Thus a careful observation of the large- 
period, and robust, Sangac interference allows, in prin- 
ciple, to extract the nanotube parameters, namely the 
interaction strength g and the velocity mismatch u from 
the Fourier transform of the conductance as a function 
of bias voltage, up to temperatures much higher than the 
Fabry-Perot oscillations temperatures. 

The temperature dependence of the Sagnac interfer- 
ence in the nanotube case is qualitatively similar to 
the single channel case. In the absence of interactions 
(g = 1), the interference can be observed to the scale T* 
proportional to vf/u; in the presence of strong interac- 
tions, however, T* becomes only weakly dependent on 
u. Unlike the single channel case, T* in the nanotube 
case is also only very weakly dependent on g in the range 
g < 0.5. This is due to the fact that only one of the four 
modes which diagonalize the Hamiltonian are interacting 
and depend on g. For the same reason, T* is higher in 
the case of the nanotube than in the single channel case, 
i.e. the reduction of T* due to interactions is not as se- 
vere in the nanotube case. The temperature dependence 
on u and g is plotted in Figure [7] In the range mentioned 
above, 23 AG « 2.8 g£«72$ P . 

IV. SUMMARY AND CONCLUSIONS 



In this paper we investigated the conductance oscilla- 
tions in carbon nanotubes due to Sagnac interference. In 
addition to theoretical interest in this large-period inter- 
ference mode, the motivation for our study also comes 
from a recent experimental realization of carbon nan- 
otube loops^. The same interference mode can arise also 
without the loop geometry in the presence of internode 
backscattering in the nanotube, as pointed out in Ref. 

m. 

The source of the Sagnac conductance oscillations is 
the difference in the velocities of left and right moving 
excitations in a carbon nanotube when the chemical po- 
tential is tuned away from half filling. Compared to the 



more familiar Fabry-Perot oscillations^, Sagnac oscilla- 
tions are expected to have a much larger period in gate 
voltage, and, as we show, in non-interacting wires sur- 
vive to a temperature a factor of vf/u higher than that 
required to observe Fabry-Perot oscillations. 

In interacting electronic wires, the above temperature 
estimation for free fermions does not apply. Our results 
for a single channel Luttinger liquid are that Tsag be- 
comes only weakly dependent on vp/u, although still 
strongly dependent on g. From our g = 0.5, 0.25 re- 
sults, the enhancement of relative to the FP interference 
is roughly: Tsag ~ 15<??Vp in the range u/vp < 0.1. 

For a strongly interacting armchair nanotube, g < 0.5, 
we find that Tsag becomes not only weakly dependent 
on vf/u, but also nearly independent of g. The Sagnac 
interference is expected to survive upto T£ AG ~ 3^jj ~ 
7 — KTpp. Considering that Fabry-Perot oscillations have 
been observed in nanotubes up to T = 10K— , Sagnac 
oscillations should be observed up to about 70K in nan- 
otubes, despite the strong interactions. 

There is also something to be learnt from examining 
the behavior of the conductance as a function of the ap- 
plied voltages. We saw that Sagnac oscillations would 
have a large period of oscillations in the applied gate 
voltage V g \ this period itself is a function of the gate 
voltage, through the dependence of the velocity differ- 
ence vr — vl = 2u. Using typical values of a nanotube 
parameters (e.g. Ref. [H), v F — 8 ■ 10 5 m/s, L = 7fim, 
g = 0.3 and a = 1/30, the period of oscillation in the gate 
voltage would be AV„ = 2 ^ F « 17V. consistent 

with the observed oscillations in Ref. [l8|. 

On the other hand, oscillations of the conductance 
as a function of the applied bias voltage V s d depend 
not only on the bare velocities, but also on the inter- 
action strength. A Fourier transform of the Sagnac os- 
cillations as a function of V s d, we show, contains four 
different frequencies corresponding to the four different 
velocities in the problem, which are roughly vp ± u,vf 
and up/ff- Using the same parameters as above we get 
AV sd = 2 ^ F w 0.5mV. This period is much smaller 
than the bandwidth of a nanotube which is a few cV, 
so in principle many oscillation periods can be observed 
and the longer period oscillations should also be mea- 
surable, allowing the slower frequency oscillations to ap- 
pear in the Fourier transform. Observation of these 
frequencies would allow us to read off the parameters 
of the nanotube, vp, u and g, at temperatures up to 
Tg AG w 70mK , which is higher than the temperatures 
associated with Fabry-Perot oscillations. 

In the single channel case, for non-interacting elec- 
trons, we were able to extract an analytic expression for 
the temperature behavior of the conductance gate volt- 
age oscillations : 

G co {T) / 2Tik B LT \ /_u \ 1 

G co (T = 0) ^ hv F J \v F J S mh(27rk B T^^) 

(32) 

and it is apparent how the ratio vp/u directly enters 



10 



the temperature scale. Unfortunately, we were so far un- 
able to extract analytic expression for T$ag m terms of g 
and u/vf for the interacting single channel or interacting 
nanotube cases, inspite of the progress on the qualitative 
understanding our numerical results allow. Such an an- 
alytical understanding should be the focus of a future 
effort. 

As can be observed in Figs. [1] and [3J the paths giving 
rise to the Sagnac intereference are similar to the paths 
that give rise to weak localization phenomena in 2d dis- 
ordered conductors. In this work we also essentially show 
that even in the presence of strong interactions, the in- 
terference survives. It is tempting to extrapolate from 
our results that weak localization should also survive 
strong interactions. This, however, is presumably true 
so long that scattering events are dominated by small 
momentum transfer. Nevertheless, our results suggest 
that a Luttinger liquid with charge and spin modes will 
still exhibit weak-localization effects, but suppressed, and 
only weakly dependent on the detuning between counter 
propogating electrons. Therefore the magnetoresistancc 
should also be strongly suppressed at low fields. 
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APPENDIX A: DIAGONALIZING THE 
HAMILTONIAN WITH NODE AND SPIN 
DEGENERACIES 

In this appendix we show how to diagonalize the 
Hamiltonian of Eq. ([24]) . where diagonalizing en- 

tails finding the appropriate change of basis that will 
transform H± c h to the sum of four Hamiltonians, each 
having a form resembling the single channel H\ c h of Eq. 
([5]). We also explain here the approximations we have 
used in our calculation. 

The first step in the diagonalization of Hb is to change 
the basis from the spin up/down to the spin symmet- 
ric/antisymmctic basis at each node: 



'a ■ = 



± 



V2 



(Al) 



applying the same transformation to the <^>'s as well. 
We notice that the density-density interaction term in- 
volves only the spin symmetric fields 9 a+ , hence the spin 
antisymmetric fields decouple and appear as two non- 
interacting (g = 1) copies of the single channel problem, 
described by the Hamiltonian iJicfc, with right moving 
velocity of Vp ± u and left moving velocities of Vp =p u. 
These are the fields labeled with j — 1 and j = 2 in Eq. 
(ESI). 



The Hamiltonian for the spin symmetric fields has a 
similar form to our starting point Hamiltonian, H^: 



2tt 



(A2) 



a=l,2 



+ (-l)° +1 2— V^ + V^ H 

Vp 



J dx2\{ ~ V9a + 



a=l, 2 



In the absence of the u term, H + is easily diagonal- 
ized by taking the node symmetric and antisymmetric 
combinations of the fields: 



± 



'3/4 



'2 + 



V2 



(A3) 



The resulting diagonal Hamiltonian would be: 



H 



+ lu=0 



' 2tt 

hv4, 

"~2T 



dx 



dx 



-(V0 3 ) 2 

93 

1 



-(W 4 ) 2 



53(V0 3 ) 2 
-<?4(V0 4 ) 2 



(A4) 



with V3 = vp, t>4 = Vp/g, <?3 = 1 and (74 = g. 

When we consider u ^ 0, it is still possible to apply 
a g and u dependent transformation to the fields, that 
will restore H + to the form in Eq. (|A4|) . with velocities 
D3/4 given by Eq. (|26p . The field mixing this transforma- 
tion entails, however, considerably complicates the book 
keeping in our perturbative calculation. Fortunately, we 
can show that a good approximation is to simply set u to 
zero in H + when the interactions are strong, and simply 
use the transformation given by Eq. (|A3|) . The first indi- 
cation that this approximation is valid is that the exact 
velocities V3/4 differ from the u = velocities by at most 
1% in the entire range of parameters we are interested 
in, which is u/vp < 0.1 and g < 0.5 

Another indication that this approximation is valid 
comes from the analysis of the single channel problem 
in Section [H] In the single channel case we derived exact 
expressions for the Sagnac interference, and found that 
for g = 0.5 and g = 0.25, the temperature dependence 
is only weakly dependent on u/vp] furthermore u only 
enters directly in the expression for the oscillation pe- 
riod of the conductance as a function of gate voltage, the 
dependence we have explicitly in our expression for the 
coherent current I co , Eq. [T5l 

Finally, wc can also calculate the exact combination 
of fields that diagonalizes H + , and verify that indeed 
they are very close to the node symmetric/antisymmetric 
combinations for the range of g and u of interest. As an 
example, the explicit change of basis from the node sym- 
metric/antisymmetric basis to the diagonalizing basis for 
g = 1/2, to second order in u/vp, is: 



144 \v F I 





H4a;4 





5 ( u 





89 / u \2 
144 I VF I 

-2^2 I u \ 
3 \v F > 





29 / it 
36 ^ v F 





22/2(_u_\ \ 
3 U F / \ 







;(^) 2 
(A5) 
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We sec that the is matrix is close to the identity matrix 
I4 X 4, since ^ C 1. The deviation from the identity 
becomes even smaller for smaller g. Note that for g m 1 
the corresponding change of basis matrix is not close to 
the identity matrix and our approximation fails. 

We stress that setting u to zero in H+ is simply a 
good numerical approximation which simplifies the cal- 
culation, and not equivalent to setting u to zero in the 
entire problem, as u still appears in spin anti-symmetric 
part of the Hamiltonian (where g = 1), and also in the 
gate voltage dependence. 

APPENDIX B: CORRELATION FUNCTIONS 

Let us now connect the explicit expressions for the co- 
herent and incoherent currents given in Section III CI and 



C e {x,t) 



R {x,t) 



where S is a short distance cutoff, i>k/l = vp/g ± u, 
and Q(x) is the step function. As mentioned in Ref. 
l26l . it is important to remember that the step functions 
are not infinitely sharp, and have a transition width of 
order a, the cutoff. The functions C^{x,t) and i?*(a;,£) 
are obtained from C e (x,t) and R e (x,t) by replacing the 
prefactor g with |. The function C e ^ is given in Eq. 

(fT2|). and: 

R H {x,t) = (B4) 
- |e(t) [Q(x)Q(x - v R t) - ®{-x)Q(\x\ - VL t)] . 

The currents are expressed in integrals over complicated 
combinations of such correlation functions. For exam- 
ple, the coherent part of the current, given by Eq. (fT5j) . 
involves the following combinations: 



C C0 (L, t) = 2C H {L, 1) - 2C e ' t '(-L, t) (B5) 
+ 2C e (0, t) - 2C«(L, 0) + C 9 (L, t) + C 9 (-L, t) 
+ 2O(L,0) - 2C*(0,t) + 0>(L,t) + C*(-L,t) 



Section IIIII Equations (JT5J) and (JT6J) , using the correla- 
tion functions defined in Section III Bl 

It is useful to define the following combination of C 6 : 
C e (x,t) = C e (0,0) -C (x,t) (Bl) 

and similarly for C^. 

In the single channel case discussed in Section [TH there 
are only a single 9 field and a single <j> field, with the 
Hamiltonian given by Eq.([5]). Since the Hamiltonian is 
quadratic we can easily evaluate all the equilibrium cor- 
relation functions at finite temperature, paying attention 
to the different time orderings that appear as a result of 
the two branches of the Keldysh contour. The results for 
finite temperature is: 

-*)]; 

(B2) 
(B3) 

I 

and 

R co (L,t)=R 84 '(L,t)-R e ' f >{-L,t) (B6) 
+ R°{Q,t)+ l -R e {L,t)+ l -R d {-L 1 t) 

-R*(0 )t ) + ±R*(L,t) + ±R*(-L,t) 

The corresponding functions for the incoherent current 
are : 

Cf nco {L, t)=± {2C e *(L, t) - 2C°+(-L, t)) (B7) 
+ 2C e (0, t) - 2C e (L, 0) + C e (L, t) + C e (~L, t) 
+ 2C*(L,0) + 2C*(0,*) - C*(L,i) - d>(-L,t) 

and 

Rinco(L,t) = R e {0,t) + ^R 8 {L,t) + ^R 8 (-L,t) 

+ R*(0,t)-±R*(L,t)-±R*(-L,t). 

In a Carbon nanotubc there are four channels, rather 
than a single one. In the non-interacting case, g = 1, all 
these channels are independent and we would recover the 
results of the single channel. Equations (|B5[) and (|B6[) 
still apply for this case. When g ^ 1, the different chan- 
nels are coupled through the interaction, and we must 
find the correct combinations of the fields 8^ and 4>i a 



i Pvl . , / t(x + v L t - iS) 

log i^r m H — w L — 



,' fiv R . f-n-(x-v R t + iS) 
l0g (^- Sm H Wn 



+ (x — > —x, t 



e(x)Q(t -—) + e(-x)ea - H) 
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which decouple and therefore diagonalizc the Hamilto- 
nian. These combinations are discussed in Appendix [A"l 
This change of basis is in general a function of u/vf and 
g, and it mixes the 9 and (f> fields, which in turn com- 
plicates the functions C co and R co further. Luckily, the 
interactions in Carbon nanotubes are strong, g « 0.3, 
and in that range, the change of basis is very close to 
the usual spin/node symmetric/antisymmetric change of 
basis. If we approximate the diagonalizing fields by these 
symmetric/antisymmetric combinations, then equations 
(|B5j) and (|B6|) would apply provided we make the follow- 
ing substitutions: 

(B8) 

3=1. .4 



Where each 0j has a different set of values for vr, vl 
and g to be used in Eq. (|B2|) . The fields B\ and 8 2 cor- 
respond to the spin asymmetric combinations, which de- 
couple from the interaction, and hence have g = 1, and 
velocities vr = vp±u and vl = vfT"- The fields 63 and 
64 both have the same left and right mover velocities, U3 
and Vi respectively, given by Eq. ((2^1) . and interaction 
parameters 1 and g, respectively. 
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